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I.  INTRODUCTION 

The  USNS  General  H.  H.  Arnold  end  USNS  General  H.  S.  Vandenberg,  operated 
by  the  Air  Force  Systems  Command,  Eastern  Space  and  Missile  Center  (ESMC) ,  are 
two  Advanced  Range  Instrumentation  Ships  (ARIS)  designed  to  gather  precision 
data  on  missile  reentry  bodies  and  penetration  aids .  The  primary  mission  of 
these  ships  Is  to  collect  metric  and  signature  data  during  the  midcourse  and 
reentry  phases  of  a  ballistic  missile's  flight.  These  instrumentation  ships 
play  a  significant  role  because  they  constitute  a  flexible  expansion  of  a 
missile  test  range  beyond  the  coverage  capabilities  of  land  based  instrumenta¬ 
tion.  Therefore,  these  ships  play  a  crucial  role  in  tracking  the  terminal 
segment  of  a  reentry  vehicle's  (RV)  trajectory  that,  very  often,  falls  in 
remote  ocean  areas. 

Although  there  are  differences,  each  ARIS  in  general,  and  H.  S.  Vandenberg 
in  particular,  is  equipped  with  advanced  instrumentation  to  gather  radar. 


telemetry,  optical,  opto-radiometric,  navigation,  and  meteorological  data. 
Typically,  during  a  trajectory  mission,  an  ARIS  collects  data  on: 

(1)  Metric  information  such  as  range,  azimuth  and  elevation  of  the  RV. 

(2)  Signature  information  such  as  radar  cross  section. 

(3)  The  ship *8  location  and  heading. 

(4)  Other  relevant  telemetry,  optical  and  meteorological  variables. 

Following  a  tracking  mission,  the  collected  data  are  processed,  off-line, 

(l.e.,  in  a  post-test  environment)  to  obtain  accurate  estimates  of  metric  and 
signature  information.  These  estimates  are  required  to  accurately  characterize 
an  RV's  flight  from  its  "penetration  point"  (l.e.,  the  point  at  which  the  RV 


iters  the  atmosphere)  to  the  "splash  point"  (l.e.,  the  point  where  the  RV 
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Engineers  and  scientists  at  the  Eastern  Space  and  Hlssile  Center  recognised 
the  need  for  improving  the  accuracy  of  metric  information  derived  from  data 
collected  by  radar  trackers  aboard  the  Instrumentation  ships.  Toward  this  end, 
they  embarked  on  a  modernization  program  aimed  at  upgrading  various  systems  on 
ARXS.  For  example,  plans  are  underway  to  upgrade  the  Timing  System,  the 
Communications  System,  the  computer  hardware,  and  so  on.  In  this  context,  the 
author  of  this  report  was  asked  TO  INVESTIGATE  THE  VARIOUS  ALTERNATIVES  THAT 
MIGHT  LEAD  TO  AN  IMPROVEMENT  IN  TRAJECTORY  TRACKING  ACCURACY  THAT  IS  SUFFICIENT 
TO  MEET  THE  MORE  STRINGENT  TEST  REQUIREMENTS  OF  CURRENT  AND  UPCOMING  MISSILE 
SYSTEMS.  This  statement  is  taken  as  the  goal  of  this  research  project. 

II.  WORK  DONE  DURING  SUMMER  1980 

During  the  Summer  of  1980,  the  principal  Investigator  spent  ten  weeks  at 
ESMC,  Patrick  Air  Force  Base,  FL  and  conducted  a  feasibility  study  of  this 
project .  The  vork  reported  here  is  a  continuation  of  the  report  submitted  to 
the  Southeastern  Center  for  Electrical  Engineering  Education.  A  copy  of  the 
above  cited  report  is  attached  herewith  for  ready  reference  as  an  appendix. 

As  per  Recommendation  (3)  of  this  report,  we  strived  to  develop  an  a-0  type 
filter  whose  performance  is  comparable  to  a  Kalman  filter.  In  the  following 
section,  we  present  a  detailed  derivation  of  the  relevant  equations. 

IK .  DERIVATION  OF  EQUATIONS 

While  tracking  a  reentry  vehicle  (RV) ,  radar  returns  in  the  form  of 
noisy  measurements  can  be  processed  to  provide  smoothed  estimates  of 
position,  velocity  and  acceleration.  Such  smoothing  can  be  accomplished 
with  a  variety  of  digital  filters  ranging  in  scope  from  a  simple  extrapolator 
to  the  complex  Kalman  filter.  The  purpose  of  this  section  is  to  derive  the 
equations  for  a  a-(J  filter  which  is  easier  to  implement  than  a  Kalman  filter 
and  yet  equivalent  in  performance  to  a  Kalman  filter  -  at  least  for  the  case 
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under  consideration. 

The  central  Idea  of  this  derivation  is  to  minimize  the  mean  square 
error  in  predicted  position  of  the  RV. 

Assumptions 

1.  For  simplicity  in  presentation,  the  coordinate  systems  and  the 
associated  transformations  are  ignored. 

2.  Only  position  and  velocity  are  Included  in  the  equations.  Acceleration 
will  be  included  in  a  subsequent  derivation. 

The  Original  Equations 
Let 

XCtj)  ■  Position  of  a  RV  at  time  t^ 

e 

X^)  “  Velocity  of  the  RV  at  time  t^ 

T  -  ti+i  -  tj,  ■  sampling  interval 

By  "sampling  interval",  we  mean  the  interval  at  which  position  and  velocity 
data  are  sampled.  We  propose  to  use  a  uniform  subscript  notation  as  follows: 

The  subscript  e  stands  for  estimated  values,  p  for  predicted  values  and  m  for 
measured  values.  Thus  Xe(tj),  for  example,  stands  for  the  estimated  position 
of  the  RV  at  time  t^.  To  make  the  meaning  of  the  adjectives,  "predicted", 
"measured",  and  "estimated"  unambiguous,  a  visualization  of  the  following 
scenario  would  be  helpful.  Typically,  during  a  tracking  mission,  we  first 
predict  the  state  (position,  velocity,  acceleration,  etc.)  of  a  RV  by  whatever 
means  (including  guesswork)  that  are  available.  Ideally,  a  complete  set  of 


trajectory  equations  in  the  Earth's  gravitational  field  would  serve  this 
purpose  well.  Since  highly  accurate  equations  of  motion  are  hard  to  come  by, 
it  is  cooaon  practice  to  base  predictions  on  past  measurements.  A  simple 
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and  widely  used  predictor  Is 


Xp  <W  -  X«  (ti>  *  T  Xe  <V 


(la) 


Xp  (tl+l5  "  Xe  (ti* 

where  the  sampling  Interval  Is  assumed  to  be  small.  Initially*  as  no 
estimated  values  are  available,  one  generally  starts  with  a  set  of  nominal 
values  for  position  and  velocity.  Then,  we  aim  the  radar  In  that  general 
vicinity  and  actually  measure  the  said  state.  Due  to  various  errors,  the 
measured  values  are  not  likely  to  be  accurate.  The  central  Idea  of  filtering 
is  to  estimate  the  state  by  using  the  predicted  and  measured  values  in  a 
suitable  filtering  formula.  For  example ,  the  present  system  on  USNS  General 
H.  S.  Vandenberg  uses  the  so-called  real-time  update  filter  (or  KTUF,  for 
short) .  A  flow  chart  of  this  filtering  process  Is  shown  In  Figure  1 .  The 
equations  characterizing  this  filter  are; 


THE  REAL  TIME  UPDATE  FILTER  (RTUF) 


IB  PPS  VICTOR  ADJ. 
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Xe(ti+lJ 


Xp(ti+1}  +  °l(Xm(ti+l)  "  Xp(ti+l)J 


(2.) 


^^i+l* 


s  <w +  4- 


”  Xp^i+1^ 


(2*) 


It  is  useful  to  note  that  if  X  *  X  or  a  -  0,  then  we  made  a  perfect 

D  P 

prediction  and  as  a  result  X  would  become  equal  to  X  .  On  the  other 

e  p 

hand,  if  a  ■  1,  the  predicted  value  is  coopletely  useless  because  the 
estimated  value  solely  relies  on  the  measured  value.  In  other  words, 
a  ■  0  corresponds  to  perfect  confidence  in  our  ability  to  predict  and 
a  ■  1  corresponds  to  perfect  confidence  in  the  accuracy  of  the  measuring 
instrument.  A  similar  interpretation  can  be  given  to  the  meaning  of  3 
in  the  second  of  the  two  equations  above.  Evidently,  a  and  3  assume  values 
that  are  somewhere  in  the  interval  (0,1).  Selection  of  appropriate  values 
for  these  parameters  is  the  crux  of  the  filtering  problem. 

In  the  present  system,  the  BXUF  calculates  the  values  of  a  and  3  by 


using  the  following  formulas, 
a  -  1  -  B2 


where 


3  -  (1  -  B)‘ 


B  .  S  B  £  B 
min 


(3a) 

(3b) 


(3c) 


where  B  .  and  B  are  arbitrarily  chosen  in  the  vicinity  of  0.9  the 
■an  “** 


acutal  value  of  B  chosen  depends  upon  the  trend  in  the  radar  data.  A  rationale 
for  this  atrategem  was  given  in  reference  [5]. 
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There  axe  several  disadvantages  to  this  approach.  First,  experience  has 
shown  that  the  RTUF  Is  not  very  effective  while  tracking  accelerating  and 
maneuvering  targets .  Even  If  the  target  is  not  actually  aaneuvering,'-  the 
uncompensated  effects  of  ship's  motion  would  make  the  target  appear  as  if  it 
is  maneuvering.  This  is  an  important  drawback  of  RTUF  for  use  on  ships. 
Finally,  errors  in  prediction  and  measurement  are  not  treated  adequately. 

The  Point  of  Departure 

In  view  of  the  above  difficulties,  we  propose  to  modify  Equation  (2) 
as  follows: 

Xe(ti+lJ  "  ^p^i+l*  +  a(ti+l)[  Xm(ti+1 
•  .  0(t. ,.)  . 

y  w  ■  yw +  -f-'vw 


Notice  that  if 

°<W  . 

a(t±)  -  a 

(5a) 

and 

^i+l* 

8^)  -  3 

(5b) 

then (Equation  (4)  reduces  to  Equation  (2).  Thus,  the  only  modification  we 

are  making  here  is  to  make  the  parameters  a  and  $  depend  on  time.  Now, 

instead  of  assigning  arbitrary  values  to  a  and  3  in  an  ad  hoc  fashion,  we 

propose  to  determine  their  values  by  explicitly  including  the  statistical  nature 

of  the  errors  in  the  algorithm.  That  is,  we  would  like  to  calculate 

«(tl+l)  and  3(ti+1)  such  that  the  expected  mean  square  error  in  Xp(tt)  is  minimized. 


-ywi 

(4a) 

-  yw 

(4b) 

sssmx  Es  aga 
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In  other  words,  \ve  would  like  to  minimize  c  [  X  (t.)]  where  e[ .]  denotes 

V  p 

the  error  in  the  bracketed  quantity  and  the  bar  Indicates  expected  value. 

This  notation  leads  to  somewhat  cumbersome  equations  and  for  expediency, 
we  henceforth  propose  to  work  primarily  with  the  position  equation  whenever 
possible,  although  both  the  equations  are  required  for  a  complete  derivation. 

Derivation  of  Sew  Equations 

With  the  «[.]  notation,  the  error  in  predicted  position  and  velocity 
are ,  from  Equation  (1) , 

«[  Xp(t±+1)]  -  c  [  Xe(t1)]  +  T  e[  Xe  (t±)  ]  +Ap(t1)  (6a) 

«[  Xp  Cti+1)l  -  «  [  Xe  (t1>3  +  Av(tj)  (6b) 

where  Ap  and  Av  are  random  variables  representing  errors  in  position  and 

velocity  and  reflect  modeling  errors,  i.e.,  errors  introduced  by  using 

Equation  (1) .  Let  us  assume  that  Ap  and  Av  have  zero  mean,  variance  of 
2  2 

0  and  O  and  a  covariance  of  u  .  (The  letters  p  and  r  in  the  preceding 
p  v  pv 

sentence  are  mnemonics  for  position  and  velocity  respectively.) 

Now  our  goal  is  to  minimize 

9  — 

«  [X  (t,.,)].  Because  X  (t.,,)  is  a  linear  combination  of  X  (t.) 
p  l+i  p  lTi  e  x 

and  X  (t.),  due  to  Equation  (1),  minimisation  of 

2  2 

t  [X  (t. .,)  ]  is  equivalent  to  the  minimization  of  e  [X.  (t.)  ]  and 
p  iri  e  x 

2  * 

<  I  X#  (tj)  ] .  So  let  us  write  the  expressions  for  these  quantities . 
Applying  the  operator  e  (.]  to  both  sides  of  Equations  (1)  and  (4)  and 


substituting  (1)  Into  (4) ,  ve  get 


€  l  Xe^i+1^  ^  ”  a  ^i+l^  e  ^  ^^i+l^  * 

+  I  1-a  (t1+1)  ]  |  c  [  Xe(tl)  ] 

+  T  6  t  Xe(t1)  ]  +  Ap(ti) 


end 


«  [  Xe(t1+1)  1  -  <  I  Xe<t±)  ]  +  Av  (ti) 


+  6(ti+l 


«  t  VW  1 


-€  [  Xe(tl)  ]  -  T  €  [  Xe(t1)  ] 


-  Ap(t±) 


} 


9. 


(7a) 


(7b) 


Squaring  and  taking  expected  values,  the  expected  mean  square  errors  are 


«2 1  i,v  ]  •  °2<t1+1>  oj  <t1+l> 


+  I  1  -  a(t1+1) 


{  «!'Vi 


)  ] 


+  2T  c  [  Xe(t1)Xe(t1)  ] 


+  T2  e2  I  Xe(t±)  ]  +  Op2(tl) 


(8a) 


and 


10. 


c2i  xe(t1+1>  i  -  ( L5W  >  <  om2  (ti+1> 

T2 


+  e2[  Xe(ti)  lj 


*‘(2B(ti+l)  1  1_0(tl+l)  1  /  T)  '  €  1  Xe<ti)Xe(ti)  1 


+[1-6  (ti+1)  ]2  2  [  xe  (tt)  ] 

-( 2  0  (ti+1)  /  T)  V  <ti> 

+  av2  (t±)  +  (02  (ci+1)  /  T2)  ap2(tl)  (8b) 

where  it  is  assumed  that  the  measurement  errors  are  uncorrelated  with 

2 

zero  mean  and  variance  of  a 

m 

To  determine  the  optimum  values  of  a(ti+1)  and  0(ti+1)  *  we  take  the 
partial  derivatives  of  Equations  (8a)  and  (8b)  with  respect  to  ct(t^+^) 
and  0  (ti+^)  and  equate  them  to  zero.  The  resulting  system  of  simultaneous 
equations  are  solved  to  get 


o(ti+l)  m  {  €  1  W  1  +  216  l  Xe(ti)  Xe(ti)  1 

«!>  >/< 


+  T2  €2  I  X,^)  ]  +  C  2  (ts 


(ti+l^ 


+  e2  [  Xe(tt)  ]  +  2Te  [  X#(t±)  X^t^  ] 


+  T2e2  [  Xe(t±)  ]  +  0p4(ti)  > 


(9a) 


and 


B(ti+1>  "  {  T  c  1  W  VV  1  +  t2«2  t  Xfe(t1)  ] 

+  H'p.ttj)  ) /l  <t1+1>  +  «2  I  *.<«!>  1 

+  2Tc  [  XeCti)  Xe(t±)  ] 

+  T2£2  [  X  (t.)  ]  +  0  2  (t.)  >  (9b) 

el  pi 

Essentially,  we  have  to  evaluate  Equations  (9a)  and  (9b)  .  Before  we  can  do 
this,  we  have  to  recursively  evaluate  e.  [  X  (t^)  Xg  (t^)  ].  Toward  this  end, 

we  multiply  Equations  (7a)  and  (7b)  (before  the  expected  values  are  taken  there) 
and  taking  the  expected  values  of  the  product  expression,  we  get 


€  1  ^‘i+l^e^i+l*  1  “  '  a(ti+l)6(ti+l)  1  T  *  °m  ^i+l5 


-  [  l-a(t1+1)  ]  [  f5(t1+1)  /  T  ]  .  Xe(ti)  ] 


+  [  1  -  a(ti+1)  1  [  1-23  (t1+1)  ]  .  c  t  Xe<ti)Xe(t1)  ] 


+  T  (  <l-a(ti+1)  ]  [  B(ti+1)  /  T  ]  .  0p  (t±) 


+  [  1  -  a(t1+1)  ]  Mpv  (tA) 


Thus,  the  current  values  of  a  and  $  are  dependent  upon 

2 

1.  The  current  variance  in  measurement  error,  0  (t. .,) 

n  iti 


i-  -'*<f.w?wfyuyin « 
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2.  Previous  estimates  of  mean  squared  error  in  smoothed  (or  estimated) 

2  ~~5  • 

position  and  velocity,  e  [  X (t.)  ]  and  c  [X  (t.)  ]  . 

c  1  6  1 

3.  Previous  estlmages  of  the  variance  and  covariance  in  predictor 
errors,  Cp2  (t±) ,  o^C^)  and  Hpy  (t±) 

IV.  PROCEDURE  FOR  IMPLEMENTING  THE  FILTER 
Assume  that  we  know 

1.  op2  (t±),  a v2  (ti)  and  ppv  (tl) 

2.  c2  I  Xfi(t1)  ]  and  e2  [  Xe(t±)  1 

Step  1.  Update  the  predictor  model  defined  in  Equations  (la)  and  (11) 

for  the  given  time  T. 

Step  2.  Estimate  the  variance  in  current  position  measurement. 

(NOTE:  No  procedure  for  this  is  given  here.  A. procedure  analogous  to  the  one 
ve  use  in  Kalman  filtering  can  be  used  here .) 

Step  3.  Obtain  estimates  of  ct(ti+j)  and  8(ti+^)  using  equations 

(9a)  and  (9b) . 

Step  4.  Obtain  smoothed  (or  estimated)  values.  X  (t...),  and  X  (t, ..) 

p  l+i  p  l+i 

from  Equations  (4a)  and  (4b)  by  using  the  current  position 
measurement . 

Step  5.  Obtain  estimates  of 

f2  1  Xe(ti+1)  1  •  ^  1  Xe  (ti+l)  J  “d 

*  *  Xe*fci+1*  Xe^i+1^  ^ 

using  the  recursive  relationships  shown  in  Equations  (8)  and  (10) 


for  use  in  next  iteration 
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Step  6.  Estimate  the  variance  and  covariance  in  predicted  values 

of  Step  1  for  use  in  next  iteration. 

(NOTE:  Same  comment  as  the  one  in  Step  2.) 

Once  the  procedure  is  started,  the  calculations  are  single  to  do. 
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IMPROVEMENT  OF  TRAJECTORY  TRACKING  ACCURACY  OF  INSTRUMENTATION  SHIPS: 

A  FEASIBILITY  STUDY 

*T 

V.  Vemuri 
ABSTRACT 

The  question  of  the  feasibility  of  Improving  metric  accuracy  of  radar 
data  obtained  from  instrumentation  ships  is  investigated.  It  is  argued  that 
major  sources  r>f  error  are  tracking,  navigation  and  stabilization.  Using  avail¬ 
able  data  as  a  guide,  it  is  argued  that  substantial  improvements  in  metric 

accuracy  are  attainable  if  the  present  auto-tracking  is  upgraded  to  on-axis 

% 

tracking  with  a  Kalman- type  filter  in  the  tracking  loop.  It  is  also  recommended 
that  a  simulation  study  be  conducted  to  gain  better  Insight  into  the  nature  of 
navigational  and  stabilization  errors.  These  two  recommendations  are  considered 
to  be  most  cost-effective  within  the  constraints  of  the  mission  under  study. 
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1.  INTRODUCTION 

The  USNS  General  H.  H.  Arnold  and  USNS  General  H.  S.  Vandenberg,  operated 
by  the  Air  Force  Systems  Command,  Eastern  Space  and  Missile  Center  (ESMC),  are 
two  Advanced  Range  Instrumentation  Ships  (ARIS)  designed  to  gather  precision 
data  on  missile  reentry  bodies  and  penetration  aids.^  The  primary  mission  of 
these  ships  is  to  collect  metric  and  signature  data  during  the  midcourse  and 
reentry  phases  of  a  ballistic  missile's  flight.  These  instrumentation  ships 
play  a  significant  role  because  they  constitute  a  flexible  expansion  of  a 
missile  test  range  beyond  the  coverage  capabilities  of  land  based  instrumenta¬ 
tion.  Therefore,  these  ships  play  a  crucial  role  in  tracking  the  terminal 
segment  of  a  reentry  vehicle's  (RV)  trajectory  that,  very  often,  falls  in 

I 

remote  ocean  areas. 

Although  there  are  differences,  each  ARIS  in  general,  and  H.  S.  Vandenberg 
in  particular,  is  equipped  with  advanced  instrumentation  to  gather  radar, 
telemetry,  optical,  opto-radiometrlc,  navigation,  and  meteorological  data. 
Typically,  during  a  trajectory  mission,  an  ARIS  collects  data  on: 

(1)  Metric  information  such  as  range,  azimuth  and  elevation  of  the  RV. 

(2)  Signature  information  such  as  radar  cross  section. 

(3)  The  ship's  location  and  heading. 

(4)  Other  relevant  telemetry,  optical  and  meteorological  variables. 

Following  a  tracking  mission,  the  collected  data  are  processed,  off-line, 

(l.e.,  in  a  post-test  environment)  to  obtain  accurate  estimates  of  metric  and 
signature  information.  These  estimates  are  required  to  accurately  characterize 
an  RV's  flight  from  its  "penetration  point"  (i.e.,  the  point  at  which  the  RV 
enters  the  atmosphere)  to  the  "splash  point"  (i.e.,  the  point  where  the  RV 
hits  the  ocean  surface). 

Engineers  and  scientists  at  the  Eastern  Space  and  Missile  Center  recognized 

the  need  for  improving  the  accuracy  of  metric  information  derived  from  data 

collected  by  radar  trackers  aboard  the  instrumentation  ships.  Toward  this  end, 

they  embarked  on  a  modernization  program  aimed  at  upgrading  various  systems  on 

2 

ARIS.  For  example,  plans  are  underway  to  upgrade  the  Timing  System  ,  the 

3 

Communications  System  ,  the  computer  hardware,  and  so  on.  In  this  context,  the 
author  of  this  report  was  asked  TO  INVESTIGATE  THE  VARIOUS  ALTERNATIVES  THAT 
MIGHT  LEAD  TO  AN  IMPROVEMENT  IN  TRAJECTORY  TRACKINC  ACCURACY  THAT  IS  SUFFICIENT 
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TO  MEET  THE  MORE  STRINGENT  TEST  REQUIREMENTS  OF  CURRENT  AND  UPCOMING  MISSILE 
SYSTEMS.  This  statement,  therefore,  Is  taken  as  the  goal  of  this  research. 

This  report,  therefore,  la  the  result  of  a  feasibility  study.  The  recom¬ 
mendations  made  at  the  end  of  this  report  constitute  a  plan  to  reach  the  goal. 

II.  OBJECTIVES 

Although  the  Instrumentation  ships  gather  a  wealth  of  data  during  a  track¬ 
ing  mission,  the  scope  of  the  goal  statement  is  confined  to  a  study  of  the 
acquisition  and  processing  of  metric  data  only;  that  is,  data  about  range, 
azimuth  and  elevation  of  a  RV.  Even  with  this  restriction,  the  scope  is  still 
too  broad  in  the  sense  that  it  allows  the  possibility  of  including  a  host  of 
competing  alternative  paths  in  reaching  the  stated  goal.  Some  of  these  alterna¬ 
tives  are  considered  extraneous  to  the  major  thrust  of  this  research  effort  and 
are,  therefore,  nominated  forthwith  from  further  consideration.  Thus,  no 
proposals  for  improving  the  electronics  aboard  the  tracking  ship  are  either  made 
or  considered,  Sims i^rly,  no  proposals  to  increase  the  number  of  tracking  ships 
are  either  made  or  considered.  In  fact  an  attempt  is  being  made  to  strike  a 
balance  between  t  desire  for  a  thorough  exploration  of  all  possibilities  and  a 
desire  to  confine  the  effort  to  those  alternatives  that  appear  feasible.  In¬ 
cluded  in  the  feasibility  considerations  are  technical  criteria  such  as  practical 
realizability  and  compatibility  with  existing  configurations,  economic  criteria 
such  as  development  and  maintenance  cost  and  the  management  criterion,  namely 
an  ability  to  meet  deadlines.  In  our  pursuit  to  reach  the  goal  within  the 
framework  of  constraints  articulated  in  the  preceding  paragraph,  we  propose  to 
seek  specific  answers  to  the  following  broadly  stated  questions: 

(1)  Can  we  improve  the  metric  accuracy  of  the  instrumentation  ships  to  meet  the 
stringent  test  requirements  of  current  and  up-coming  missile  systems? 

(2)  If  the  answer  to  this  question  is  in  the  affirmative,  then  to  what  extent 
can  the  accuracy  be  improved  with  the  present  configuration  of  the  ship's 
instrumentation? 

(3)  If  modifications  to  the  present  configuration  are  warranted,  then  can  these 
modifications  be  made  within  reasonable  time  and  cost  constraints? 

(A)  If  modifications  are  considered  necessary,  which  of  these  modifications 
should  be  directed  at  improving  the  quality  of  the  data  acquisition  phase 
and  which  toward  improving  the  quality  of  the  post-test  data  analysis  phase? 
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(5)  Whatever  the  modifications  are,  what  would  it  take  to  actually  implement 

them? 

III.  DESCRIPTION  OF  A  TYPICAL  ARIS  MISSION 

To  present  the  results  of  this  research  in  a  proper  perspective,  it  is 
useful  to  review  the  salient  steps  of  a  typical  mission. 

Stated  in  a  nutshell,  the  present  tracking  procedure  consists  of  two  broad 
phases:  the  on-board  data  acquisition  phase  and  the  off-board  data  processing 
phase.  Preparatory  to  the  data  acquisition  operations,  first  the  longitudinal 
axis  of  the  ins trumentation  ship  is  positioned  roughly  perpendicular  to  the 
plane  of  the  trajectory  of  the  RV.  Then  the  location  of  the  ship  with  respect 
to  an  earth-fixed  geocentric  (EFG  coordinate)  system  is  determined  with  the  help 
of  data  provided  by  the  ship's  inertial  navigation  system  (SINS)  as  well  as  from 
fixes -obtained  from  an  array  of  submerged  transponders.  Then,  to  lend  a  degree 
of  stability  to  the  observation  platform,  the  ship  is  maintained  on  a  linear 
course  at  a  small  constant  velocity.  Then  the  target  is  acquired  by  the  radar 
and  tracked  in  the  so-called  "auto-track"  mode  and  the  necessary  data  are 
recorded  on  digital  and  video  tapes.  This  constitutes  the  end  of  the  first  phase. 

Among  all  the  data  items  collected  during  the  data  acquisition  (or  first) 
phase,  only  a  few  are  of  particular  interest  to  us  here.  Those  data  items  are 
navigation  data  (to  determine  the  exact  location  and  attitude  of  the  ship  at 
any  time),  metric  data  (to  determine  the  range,  elevation  and  azimuth  of  the 
target  with  respect  to  any  desired  coordinate  system)  and  timing  data  (for 
synchronization  purposes).  These  data  items,  recorded  on  digital  tapes,  will 
be  referred  to  as  raw  mission  data. 

The  second  phase  of  the  tracking  procedure,  called  the  data  processing  phase, 
takes  the  raw  mission  data  as  input  and  produces,  as  output,  target  position, 
velocity  and  acceleration  relative  to  any  desired  fixed  point  and  any  fixed 
reference  system.  That  is,  the  output  of  this  phase  is  the  required  metric  data 
about  the  RV's  trajectory.  This  phase  is  conducted  off-line  in  a  post-test 
environment. 

IV.  DISCUSSION  ON  PRESENT  EFFORTS  TO  IMPROVE  METRIC  ACCURACY 

The  purpose  of  this  section  is  to  identify  and  coa&ent  on  sosw  of  the  problem 
being  encountered  in  the  context  of  accuracy  improvement.  Engineers  at  ESMC 


(representing  USAF,  Pan  Am,  RCA)  are  addressing  these  problems  systematically 
and  thoroughly.  Their  approach  is  based  on  an  identification  of  all  possible 
sources  of  error,  classification  into  categories  (significant,  insignificant; 
biased,  random;  and  so  on),  and  charaterization  in  terms  of  variables  that  can 
be  controlled.  This  approach  led  to  a  satisfactory  charaterization  of  the 
behavior  of  some  of  the  errors.  For  example,  engineers  at  ESMC  appear  to  be 
fairly  happy  with  their  ability  to  control  the  errors  due  to,  say,  nonortho¬ 
gonality  of  antenna  axes,  antenna  droop,  refraction  and  RF  optical  errors. 

This  is  understandable  because  such  errors  are  common  with  all  (land  based  and 
shipboard)  radar  trackers  and  most  of  the  technical  personnel  have  rich  experi¬ 
ence  in  this  area. 

I 

The  problem  gets  a  little  complex  while  studying  the  performance  of  ship¬ 
board  radars.  A  new  class  of  errors  enter  the  picture  here.  A  case  in  point 
are  errors  introduced  due  to  our  inability  to  estimate  accurately  the  position 
add  attitude  of  the  ship  on  the  ocean  surface.  Position  simply  means  the 
coordinates  of  the  ship  in  the  EFG  (earth  fixed  geocentric)  coordinate  system. 
Attitude  refers  to  two  aspects  of  ship  motion:  roll  and  pitch  relative  to  the 
local  vertical,  and  heading  relative  to  north.  In  this  report,  we  would  like 
to  refer  to  position  and  heading  errors  as  navigational  errors  and  errors  due 
to  roll  and  pitch  as  stabilization  errors.  In  addition  to  these  errors,  we 
also  have  to  deal  with  flexure  errors,  that  is,  errors  caused  by  bending  and 
flexing  of  the  ship's  body.  Some  effort  is  already  under  way  to  minimize  the 
effect  of  these  errors.  For  example,  attempts  are  being  made  to  minimize 
navigational  errors  by  supplementing  SINS  with  position  data  obtained  from 
submerged  transponders.  Similarly,  attempts  are  being  made  to  minimize  flex¬ 
ure  errors  by  redesigning  the  physical  layout  of  some  of  the  instruments. 

This  Investigator  feels  that  all  the  above  cited  measures  are  necessary 
but  not  sufficient.  This  feeling  is  based  on  the  assumption  that  there  always 
exists  residual  errors;  to  assume 'that  we  can  identify,  characterize,  and  com¬ 
pensate  each  and  every  source  of  residual  error  is  unrealistic.  There  is  no 
suggestion  here  to  indicate  that  the  present  approach  of  correcting  major 
errors  on  an  individual  basis  be  abandoned.  What  is  being  suggested  here  is 
to  supplement  the  present  effort  with  a  good  tracking  procedure  that  explicitly 
recognizes  the  need  to  compensate  the  residual  errors  during  the  data  acquisi¬ 
tion  phase. 
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the  development  of  good  tracking  procedures  for  shipboard  radars  Is  not 
a  trivial  matter.  The  difficulty  arises  because  the  errors  associated  with 
the  range,  azimuth  and  elevation  (say  in  the  EFG  coordinate  system)  are  an 
aggregate  of  residual  errors  from  all  sources  cited  above.  To  gain  an 
insight  into  the  nature  and  magnitude  of  these  aggregate  residual  errors. 

It  Is  useful  to  digress  and  inspect  the  magnitudes  of  these  errors  resulting 
from  the  tracking  procedures  currently  in  use  at  ESMC. 

The  present  data  acquisition  system  on  H.  S.  Vandenberg,  for  example, 
uses  the  autotracking  procedure.  The  monopulse  radar  in  an  autotrack  mode 
is  essentially  a  simple  closed  loop  control  system  in  which  the  radar  receiver 
feeds  the  control  signal  directly  to  the  radar  sensors.  In  other  words,  there 
is  no  predictive  capability  in  the  autotracking  mode.  Improved  performance 
can  be  obtained  if  the  loop  "is  ‘closed  through  a  computer  so  that  the  computer 
can  be  used  to  anticipate  the  RV*s  position  from  past  data.  This  is  the  basic 
idea  of  on-axis  tracking.  The  Eastern  Test  Range  used  this  on-axis  concept  as 
early  as  1967  to  control  the  ground  based  radar  (the  AN/FPQ-13)  at  Grand  Bahama 
Island.  The  Western  Test  Range  also  used  the  on-axis  concept  for  its  ground 
based  radar  at  Kaena  Point,  Hawaii.4 

The  on-axis  concept  can  be  used  at  several  levels  of  sophistication.  Even 
the  most  primitive  type  of  on-axis  tracking  is  an  order  of  magnitude  better 
than  the  simple  autotracking.  Indeed,  the  on-axis  procedure  used  in  Grand 
Bahama  Island  and  Kaena  Point  radar  trackers  (namely,  the  *-8  tracking)  is 
relatively  primitive.  Nevertheless,  it  is  apparently  giving  satisfactory  per¬ 
formance.  For  example,  it  is  claimed  that  the  above  two,  as  well  as  other  ground 
based  radars  using  *-8  tracking,  are  giving  an  accuracy  of  +  6  ft.  in  linear 
measurements  and  +  (1/2  )  x  2tt  radians  (=  0.006  milliradians)  in  angle  measure¬ 

ments.^  In  contrast,  the  autotracked  radars  on  H.  S.  Vandenberg  are  committing 
errors  of  up  to  +  1500  ft.  in  range  and  +  0.5  milliradians  in  angle  measurements.^ 
Two  possible  reasons  for  the  poor  quality  of  metric  information  from  H.  S. 
Vandenberg  are:  (a)  measurements  are  taken  from  8  moving,  unstable  platform, 
namely  the  ship;  and  (b)  the  radar  trackers  on  the  ship  are  using  autotracking 
rather  than  on-axis  tracking.  A  useful  question  that  should  be  answered  at 
this  stage  is:  What  part  of  the  error  is  caused  by  the  moving,  unstable  plat¬ 
form  and  what  part  by  auto tracking?  For  example,  1  mr  of  ship  roll,  it  is 
observed,  produces  about  one  inch  of  C-band  antenna  movement  but  can  introduce 
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up  to  1000  ft  of  RV  position  error  at  pierce  point. ^  However,  it  is  dangerous 
to  conclude  (from  what  is  stated  so  far)  that  2/3  of  the  range  error  is  due  to 
an  unstable  platform  and  only  1/3  due  to  autotracking;  other  factors,  such  as 
-navigational  errors,  also  play  a  role  in  this  context.  For  example,  it  appears 
that  the  uncertainty  in  this  ship's  location  could  be  of  the  order  of  +  1500  ft. 
The  precise  impact  of  this  error  on  metric  information  is  unclear. 

This  digression  (the  preceding  three  paragraphs)  confirms  the  general 
qualitative  feeling  that  tracking  stabilization  and  navigational  problems  are 
major  problems  to  be  solved.  However,  ve  do  not  yet  have  any  quantitative 
measures  of  the  contribution  of  each  of  these  to  the  total  observed  error. 

Until  we  can  develop  quantitative  estimates  of  the  error  contributions  from 
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each  of  these  sources,  we  will  not  be  in  a  strong  position  to  decide  where  we 
should  put  our  resources  to  meet  our  goal.  Toward  this  end,  we  developed  in 
the  next  section,  some  quantitative  estimates  of  the  order  of  magnitude  improve¬ 
ments  that  can  be  made  by  improving  tracking.  We  do  not  have  information  to 
help  us  make  analogous  estimates  insofar  as  navigation  and  stabilization  errors 
are  concerned.  In  Section  VI,  we  proposed  that  a  simulation  study  be  conducted 
to  gain  better  Insight  into  navigation  and  stabilization  errors. 

V.  SPECTRUM  OF  FUTURE  POSSIBILITIES 

The  limited  scope  of  the  discussion  in  the  preceding  four  sections  is  not 
sufficient  to  completely  reveal  the  complexity  of  a  tracklngmission.  Neverthe¬ 
less,  the  success  of  a  mission  depends  upon  the  coordinated  operation  of  a 
number  of  subsystems.  Modification  of  one  subsystem  or  one  operational  pro¬ 
cedure  could  create  a  ripple  effect  demanding  modifications  in  other  systems. 

As  the  scope,  as  outlined  in  Section  II,  of  this  effort  is  somewhat  limited,  it 
is  not  possible  to  address  all  these  ripple  effects.  By  the  same  token,  they 
cannot  be  ignored  either.  In  this  section,  we  intend  to  propose  some  modifica¬ 
tions  and  briefly  touch  upon  some'  of  the  possible  ripple  effects. 

A.  Possible  Modifications  to  Tracking.  The  radar  trackers  on  H.  S. 

Vandenberg  now  use  autotracking.  Engineers  at  ESMC  recognized  a  decade  ago 

8  9 

the  inadequacy  of  autotracking  and  the  desirability  of  on-axis  tracking.  * 

In  fact.  Reference  9  gives  the  figures  shown  in  Table  1  to  illustrate  the 
advantage  of  on-axis  tracking  over  autotracking.  The  reference  did  not  mention 
the  filter  used  in  the  on-axis  algorithm.  However,  there  is  reason  to  believe 
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that  it  probably  was  an  «-8  filter;  definitely  not  a  Kalman  filter. 

Table  1 

A  comparison  of  orbital  residuals  from  echo  track  of  the  satellite 
Pegasus  performed  by  Radar  12.15  at  Ascension  in  February  1974. 


1 

Absolute 

Z  Improvement 

Variable 

Autotracking 

On-Axis 

Improvement 

(Approximate) 

Elevation  errors 

0.22  mr 

0.044  mr 

0.176 

802 

Range  error 

33.5  ft 

10.5  ft 

23.0 

692 

inspection  of  Table  1  reveals  that  accuracy  improvements  of  the  order  of 
70,  -  80Z  can  be  attained  by  replacing  auto tracker  with  a  simple  «-0  tracker. 
Furthermore,  in  Reference  4  (page  1.5),  it  was  estimated  that  a  nine-variable 
Kalman  filter  has  the  potential  of  giving  a  302  improvement  in  real-time 
tracking  accuracy  over  an  «-8  filter.  Admittedly,  these  two  estimates  were 
made  in  the  context  of  ground  based  radars.  In  spite  of  this  knowledges,  it 
is  surprising  that  ESMC  did  not  initiate  plans  to  upgrade  tracking  procedures 
on  ARIS  until  the  1980s.  There  are  two  frequently  quoted  reasons  for  this 
delay.  First,  that  the  ARIS  were  originally  designed  to  gather  signature  data 
and  the  desire  to  use  them  for  metric  data  gathering  is  an  afterthought. 

Second,  that  everything  possible  is  being  done  to  correct  errors  due  to  various 
sources  and  that  there  is  nothing  much  to  be  gained  by  overhauling  the  tracking 
procedure.  The  first  of  the  above  implies  that  there  is  a  change  in  mission 
requirements  and  the  second  reflects  a  lack  of  strong  conviction.  Nevertheless, 
plans  are  under  way  to  explore  the  possibility  of  replacing  autotracking  with 
on-axis  tracking.  This  is  definitely  a  step  in  the  right  direction, 
although  not  a  decisive  one.  The  author  feels  that  this  line  of  thinking  should 
be  pursued  more  aggressively  by  going  all  the  way  to  a  Kalman-type  tracker  rather 
than  stop  in  midstep  with  an  «-0  tracker.  If  there  are  some  practical  reasons 
for  a  reluctance  to  go  all  the  way  for  a  full-blown  Kalman  tracker,  there  exist 
some  suboptlmal  Kalman  trackers  for  consideration.  There  are  even  optimal  *-8 
trackers  that  are  equivalent  to  Kalman  trackers  in  their  performance.*® 
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Without  going  into  detailed  analyses  of  the  advantages  and  disadvantages 
of  Kalman  filters  in  general,  let  us  briefly  look  into  the  ripple  effects  o*f 
on-axis  traking  in  the  context  of  a  typical  ARTS  mission.  At  present,  most 
of  the  metric  information  is  being  derived  by  off-line  processing  of  data  in 
a  post  test  environment.  On-axis  tracking  (be  it  «-($  or  Kalman  variety)  is  an 
on-line,  real-time  procedure.  This  implies  that  there  is  a  need  for  an  addi¬ 
tional  computer,  on-board,  dedicated  to  on-axis  tracking.  This  computer  should 
be  supplied  with  navigation  data,  in  addition  to  the  usual  tracking  data.  This, 
in  turn,  means  that  we  cannot  afford  the  luxury  of  waiting  for  the  splash  point 
to  occur  in  order  to  determine  the  ship's  location.  That  is,  we  have  to  make 
navigation  independent  of  tracking  but  not  vice  versa. 

A  second  possible  side  effect  has  to  do  with  data  processing  operations 
such  as  the  editing  of  raw  data  and  accuracy  of  encoders.  For  example,  Kalman 
filters  tend  to  be  sensitive  to  the  editing  scheme  used.  Also,  if  a  lot  of 
effort  is  going  to  be  expended  in  improving  the  accuracy  of  data  gathered, 
corresponding  attention  should  be  paid  in  maintaining  this  accuracy  in  digitiz¬ 
ing  and  encoding  this  information.  Toward  this  end  the  UNIVAC  1219  type  computer 
with  a  standard  word  length  may  not  be  sufficient.  A  computer  with  a  32-bit  word 
length  and  64K  of  storage  is  probably  needed  for  any  sophisticated  scheme. 

A  third  side  effect  of  on-axis  tracking  is  the  need  to  perform  the  computa¬ 
tions  on-line  although  this  is  not  at  all  a  requirement  of  the  mission. 

B.  Possible  Modifications  to  Navigation.  There  is  a  general  feeling  at 
ESHC  that-  the  impact  of  navigational  errors  on  metric  accuracy  art  «->t  as  severe 
as  those  of  tracking  and  stabilization  errors.  This  is  predicated  on  the 
assumption  that  the  coordinates  of  the  splash  point  can  always  be  determined 
fairly  accurately  and  the  ship's  location  with  respect  to  the  splash  point  can 
therefore  be  derived  because  the  ship  is  generally  not  too  far  from  the  splash 
point.  As  noted  earlier,  this  implies  that  the  ship's  location  can  be  accurately 
determined  after  the  mission,  not  before  or  during  the  mission  because  the  splash 
point  is  the  terminal  point  of  the  RV  trajectory.  But  the  accuracy  of  on-axis 
tracking  depends,  to  some  extent,  on  advance  knowledge  of  ship's  position.  Thus, 
to  make  the  on-axis  tracking  really  useful,  we  must  look  for  ways  to  determine 
the  ship *8  position  and  heading  by  methods  that  do  not  depend  upon  data  derived 
from  the  trajectory.  This  is  an  Interesting  side  effect  of  using  on-axis  track¬ 
ing.  Therefore,  either  we  have  to  determine  ship's  position  only  from  ship's 
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Inertial  navigation  system  (SINS)  or  supplement  it  with  some  other  system.  One 

possibility  is  to  explore  the  possibility  of  using  the  Navstar  Global  Positioning 

System  (GPS)  currently  under  development  by  the  Department  of  Defense.  The  GPS 

system  can  be  used  to  get  both  time  and  position  information  precisely.  In  fact, 

2 

the  GPS  system  is  being  considered  in  Timing  Modernization  Plan.  If  it  is  going 
to  be  used  to  get  timing  information,  one  may  as  well  use  it  to  get  position  data 
also.  Presumably,  the  cost  involved  in  using  the  GPS  system  is  in  building  a  GPS 
receiver. 


Alternatively,  one  can  take  advantage  of  the  similarity  between  the  problem 

of  tracking  the  position  of  an  RV  in  space  and  the  problem  of  tracking  the  position 

of  a  ship  on  ocean  surface.  As  both  problems  involve  the  observation  of  a  moving 

object  with  imperfect  instruments  and  subsequently  filtering  the  noisy  data,  the 

12-14 

same  procedure,  with  appropriate  modification,  can  be  used  in  both  cases. 

This  strategy  has  the  aesthetic  appeal  of  depending  on  a  uniform  procedure  to 
splve  a  broad  class  of  problems  in  the  mission.  Such  a  streamlined  procedure, 
that  keeps  the  number  of  new  things  to  be  learned  to  a  minimum,  vastly  improves 
the  efficiency  of  people  who  design  and  maintain  these  facilities. 

C.  The  Platform  Stabilization  Problem.  A  possible  solution  to  the  platform 
stabilization  problem  appears  to  be  a  little  more  difficult.  Although  Kalman-type 
approaches  were  proposed  in  the  past  by  several  investigators  to  solve  analogous 
problems,  their  applicability  to  the  present  problem  needs  to  be  Investigated.  A 
simulation  approach  appears  to  offer  an  ideal  compromise  here.  The  simulator 


would  take  sea  state,  wind  and  ship  velocity  as  inputs  and  produces  expected  pitch, 
■  ,  ‘  15»  16 

roll  and  yaw  as  outputs.  The  results  of  this  simulation  can  be  used  to 

predict  errors  due  to  ship  motion. 

VI.  RECOMMENDATIONS 

Analyses  in  the  preceding  sections  can  now  be  used  to  answer  the  questions 
raised  In  Section  II  of  this  report. 

(1)  It  is  possible  to  Improve  the  accuracy  of  metric  measurements  by  a  substan¬ 
tial  margin.  There  is  sufficient  evidence  in  published  literature  for  this 
possibility. 

(2)  If  the  objective  is  to  achieve  as  much  of  this  Improvement  as  possible  with 
minimum  change  in  the  ship's  Instrumentation,  then  the  best  course  of  action 
is  to  replace  auto tracking  with  on-axis  tracking.  This  option  would 
probably  require  a  minicomputer  dedicated  to  tracking.  Although  an  exact 
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estimate  of  the  size  of  the  computer  depends  upon  the  type  of  filter  used 
in  the  on-axls  tracking,  an  educated  guess  would  be  a  machine  with  48-64K 
of  memory  size  with  32-bit  word  length. 

(3)  If  the  on-axis  tracking  idea  is  acceptable,  then  this  investigator  feels 
that  a  serious  effort  at  implementing  a  Kalman-type  filter  be  initiated. 

In  trying  to  implement  a  Kalman-type  filter,  consideration  should  be  given 
to  a  determination  of  the  type  of  suboptimal  filter  that  best  suits  the 
needs  of  this  project.  In  view  of  the  time  and  cost  constraints,  the 
possibility  of  developing  an  optimized  -  filter  (that  is,  one  whose  per¬ 
formance  is  equivalent  to  a  Kalman  filter)  for  possible  implementation  on 
modern  microcomputers  should  not  be  Ignored. 

t 

(4)  Once  the  above  idea  is  adopted,  it  is  important  to  realize  that  now  we  are 
dealing  with  the  so-called  embedded  computer  systems.  The  reliability  and 
maintainability  of  such  embedded  computer  systems  very  much  depends  upon 

1  the  quality  of  software  support.  Toward  this  end,  it  is  strongly  recom¬ 
mended  that  all  future  software  development  efforts  follow  modern  ideas 
of  software  engineering;  that  is,  ideas  such  as  top-down  design  and 
structured  programming. 

(5)  The  author  also  feels  that  the  short  range  objective  of  improving  the 
quality  of  metric  data  from  R.  S.  Vandenberg  should  not  be  allowed  to 
cloud  the  long  range  objective  of  maintaining  the  concept  of  using  ARIS 
as  a  flexible  expansion  of  a  missile  test  range.  In  this  context,  it  is 
extremely  useful  to  conduct  a  simulation  study  to  determine  how  the  various 
aspects  of  a  ship's  motion  influence  the  overall  accuracy. 

(6)  Finally,  regardless  of  whatever  action  is  being  taken  insofar  as  ARIS  are 
concerned,  it  is  important  that  all  personnel  concerned  with  RV  tracking 
at  ESMC  be  encouraged  to  get  abreast  with  the  developments  in  technology. 
Recent  advancements,  and  some  not  so  recent,  in  guidance  and  control, 
computer  simulation,  distributed  processing  and  software  engineering  are 
revolutionizing  thinking  in  these  areas. 
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